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Abstract 

We have searched for lepton flavor violating r decays with a pseudoscalar meson (n, n' and tt°) 
using a data sample of 401 fb _1 collected with the Belle detector at the KEKB asymmetric-energy 
e + e~ collider. No evidence for these decays is found and we set the following upper limits on the 
branching fractions: £>(r~ — » e~n) < 9.2 x 10~ 8 , B(t~ — » fj,~rj) < 6.5 x 1CT 8 , £>(t~ — » e~n') < 
1.6xl(T 7 , B(t~ -> p-rf) < 1.3X10 -7 B(t~ -> e'n ) < 8.0xl0~ 8 and£(r~ -> /^vr ) < 1.2xl0- 7 
at the 90% confidence level, respectively. These results improve the previously published limits by 
factors from 2.3 to 6.4. 

PACS numbers: ll.30.Fs; 13.35.Dx; 14.60.Fg 
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INTRODUCTION 



Lepton flavor violation (LFV) is allowed in many extensions of the Standard Model 
(SM), such as Supersymmetry (SUSY) and leptoquark models. In particular, lepton flavor 
violating decays with a pseudoscalar meson (M° = 77, 7/ and 7r°) are discussed in models with 
Higgs-mediated LFV processes pj, heavy singlet Dirac neutrinos 0, R— parity violation in 
SUSY S 0], dimension-six effective fermionic operators that induce t — fi mixing and 
others p. The best upper limits for these modes are in the range (1.5—10) x 10~ 7 at 
the 90% confidence level. These were obtained by the Belle experiment using 154 fb _1 of 
data 0]. 

In this paper, we report a search for lepton flavor violating decays with a pseudoscalar 
meson r~ — > £~M° (£ = e or \i and M° = 77,77' or n° )[f] using 401 fb _1 of data collected 
at the T(4S) resonance and 60 MeV below it with the Belle detector at the KEKB e + e~ 
asymmetric-energy collider (if. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL), all located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K® mesons and to identify muons 
(KLM). The detector is described in detail elsewhere 0. 

Particle identification is very important in this measurement. We use particle identifi- 
cation likelihood variables based on the ratio of the energy deposited in the ECL to the 
momentum measured in the SVD and CDC, the shower shape in the ECL, the particle 
range in the KLM, the hit information from the ACC, the measured dE/dX in the CDC 
and the particle time-of-fli ght from the TOF. For lepton identification, we form likelihood 
ratios V(e) [l(| and V(fi) fill ] based on the electron and muon probabilities, respectively, 
which are determined by the responses of the appropriate subdetectors. 

In order to determine the event selection requirements, we use the Monte Carlo (MC) 
samples. The following MC programs have been used to generate background events: KO- 



RALB/TAUOLA 



3 for r+r", QQ for BB and continuum, BHLUMI Q for Bhabha 
events, KKMC [15] for e + e~ — > and AAFH 0] for two-photon processes. Signal MC 

is generated by KORALB/TAUOLA. Signal r decays are two-body and assumed to have a 
uniform angular distribution in the r lepton's rest frame. All kinematic variables are calcu- 
lated in the laboratory frame unless otherwise specified. In particular, variables calculated 
in the e + e~ center-of-mass (CM) frame are indicated by the superscript "CM". 



EVENT SELECTION 

We search for t + t~ events in which one r decays into a lepton and a pseudoscalar meson 
on the signal side, while the other r decays into one charged track with a sign opposite to 
that of the signal-side lepton and any number of additional photons and neutrinos on the 
tag side. Thus, the experimental signature is: 



[|] Unless otherwise stated, charge conjugate decays are included throughout this paper. 
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{r~ -> e~ or fi~) + M°(= 77,77' or 7r )} + |r + -> (a track) + + (n 7 AG > 0) + X ( missing) j. 

We reconstruct a pseudoscalar meson in the following modes: 77 — > 77 and 7r + 7r~7r (^ 77), 
77' — > p(— > 7r + 7r~)7 and r](— > 77)7r + 7r~, 7r° — > 77. While the 7r° — > 77 and 77 — > 77 modes 
correspond to the 1-1 prong configuration, the other modes give 3-1 prong configurations. 
All charged tracks and photons are required to be reconstructed within the fiducial volume, 
defined by —0.866 < cos# < 0.956, where 9 is the polar angle with respect to the direction 
opposite to the e + beam. We select charged tracks with momenta transverse to the e + beam 
Pt > 0.1 GeV/c while the photon energies must satisfy the requirement > 0.1 GeV (0.05 
GeV) for the 1-1 prong (the 3-1 prong) configuration. 

Candidate r-pair events are required to have two and four tracks with a zero net charge, 
for the 1-1 and 3-1 prong configurations, respectively. Event particles are separated into 
two hemispheres referred to as the signal and tag sides using the plane perpendicular to the 
thrust axis |17| . Whereas the tag side contains a single track, the signal side contains one 
or three tracks. For the 1-1 prong configuration, we require that the number of photons 
on the signal side be two or three. The track on the signal side is required to satisfy the 
lepton identification selection. The electron and muon identification criteria are V(e) > 0.9 
with p > 0.7 GeV/c and V(fi) > 0.9 with p > 0.7 GeV/c, respectively. The efficiencies for 
electron and muon identification after these requirements are 92% and 88%, repectively. To 
reduce fake pseudoscalar meson candidates, we reject radiative photons from electrons on 
the signal side if cos# e7 > 0.99. 

To ensure that the missing particles are neutrinos rather than photons or charged particles 
that fall outside the detector acceptance, we impose additional requirements on the missing 
momentum vector, p m i SS , calculated by subtracting the vector sum of the momenta of all 
tracks and photons from the sum of the e + and e~ beam momenta. We require that the 
magnitude of p m i SS be greater than 0.4 GeV/c and that its direction point into the fiducial 
volume of the detector. Since neutrinos are emitted only on the tag side, the direction of 
Pmigg should lie within the tag side of the event. The cosine of the opening angle between 
p m i SS and the thrust axis (on the signal side) in the CM system, cos6^ a g_ thrust , is therefore 
required to be less than —0.55. 

Event selection for the ?] — > 77 mode 

The 77 meson is reconstructed by combining two photons. The mass window is chosen 
to be 0.515 GeV/c 2 < m 77 < 0.570 GeV/c 2 , which corresponds to —3.0 and +2.5 standard 
deviations (a) in terms of the mass resolution. To avoid fake 77 candidates, we reject those 
photons that form 7r° candidates in association with any other photon with E 1 > 0.05 GeV, 
within the 7r° mass window, 0.10 GeV/c 2 < M 77 < 0.16 GeV/c 2 . To suppress fake 77 events 
from beam background and initial state radiation (ISR), we require that the higher and 
lower energy photons in an 77 candidate (E 7 i and E^) satisfy the requirement E^\ > 0.6 
GeV and E^ > 0.25 GeV, respectively. To reduce background from Bhabha and fi + fi~ 
events, we require the momentum of a lepton and a tag-side charged particle to be less than 
4.5 GeV/c. 

The total visible energy in the CM frame, E™, is defined as the sum of the energies of 
the 77 candidate, the lepton, the tag-side track (with a pion mass hypothesis) and all photon 
candidates. We require E™ to satisfy the condition: 5.29 GeV < E™ < 10.0 GeV. To 
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reduce background from two-photon and Bhabha events, we further require E™ to 

satisfy the veto condition: E™ > 8.5 GeV for the muon mode (electron mode) if the track 
on the tag side is a muon (electron). The cosine of the opening angle between the lepton and 
the rj in the CM system, cos9™ , is required to lie in the range 0.50 < cos < 0.85. The 
reconstructed mass on the tag side using a track (with a pion mass hypothesis) and photons, 
m tag , is less than 1.777 GeV/c 2 . In order to suppress background from qq (q = u,d,s,c) 
continuum events, the following requirement on the number of photon candidates on the tag 
side is imposed: n TAG < 2. 

The correlation between the momentum of the track on the tag side, p™, and the cosine of 
the opening angle between the thrust and missing particle, cos ^thnist-miss ^ n ^ ne CM system is 
employed to further suppress backgrounds from generic t + t - and p + p~ events via the follow- 
ing requirements: pgf > 1.1 log(cos0 c ^ t _ miss + 0.92) + 5.5, andp™ < 5cos^ ust _ miss + 7.8 
where p™ is in GeV/c (see Fig. Finally, we require the following relation between the 
missing momentum p m i SS and missing mass squared m^ iss to further suppress background 
from generic t + t~ and continuum background. In signal events, two neutrinos are included 
if the r decay on the tag side is leptonic decay, while one neutrino is included if the r decay 
on the tag side is a hadronic decay. Therefore, we separate events into two classes according 
to the track on the tag side: leptonic or hadronic. We apply the following requirements 
Pmiss > — 10 m miss + 4 an d Pmiss > l-l m miss — 0-3 f° r a leptonic mode on the tag side, and 
require p m i SS > — 5m^ iss — 0.25 and p m i SS > 2.1m 2 niss — 0.3 for a hadronic mode on the tag 
side, where p tag is in GeV/c (see Fig. EJ). Following all the selection criteria, the signal 
detection efficiencies for the r~ — > e~77(— > 77) and t~ — > /j,~r](-^ 77) modes are 5.08% and 
7.13%, respectively. 

Event selection for the ?] — > tt + tt~tt° mode 

The rj meson is reconstructed from 7r + 7r~7r°. The 7r° candidates are formed from a pair of 
photons that satisfy 0.115 GeV/c 2 < M 77 < 0.152 GeV/c 2 (±2. 5a). with p n o > 0.1 GeV/c, 
where p n o is the ir° momentum in the laboratory system. The mass window for rj — > 7r + 7r~7r° 
is chosen as 0.532 GeV/c 2 < M 3n < 0.562 GeV/c 2 , which corresponds to ±3. Oct. Figure El 
shows the mass distribution for the rj — > 7r + 7T~7r° candidates. Good agreement between data 
and the MC expectation is observed. 

Similarly for the rj(— > 77) mode, we require the following: 5.29 GeV < E™ <10.0 GeV, 
0.50 < cos0^ < 0.85 (see Fig. H (a) and (b)) and m tag < 1.777 GeV/c 2 . The requirement 
on the number of the photon candidates for the signal is n SIG < 1. In addition we apply 
the following requirements: p m i SS > — 5m 2 liss — 0.25 and p m i SS > 2.1m 2 nisg — 0.3 for hadronic 
tags, and p m i SS > — lOra^iss — 4 and p m i SS > l.lm^ — 1 for leptonic tags, respectively. After 
all the selection criteria, the signal detection efficiencies for the r~ — > e~ri(—> 7r + 7r _ 7r°) and 
t~ — > fJi~r](— > 7r + 7r~7r ) modes are 5.25% and 7.60%, respectively. 

Event selection for the 7/ — > p(— > tt + tt^)^/ mode 

For the p — > tt + 7t~ selection, the mass window is chosen to be 0.550 GeV/c 2 < < 
0.900 GeV/c 2 . We reconstruct 77' candidates using a p candidate and a photon on the signal 
side. The 77' mass window is chosen to be 0.930 GeV/c 2 < m pl < 0.970 GeV/c 2 , which 
corresponds to —3.0 and +2.5<r. Furthermore, we veto photons from 7r° candidates in order 
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(c) generic rr (d) fifi 




-1 -0.8 -0.6 -1 -0.8 -0.6 

COS 1 miss-thrust COS ^miss-thrust 



FIG. 1: Scatter-plots of (a) signal MC (r — > \i r}{— > 77)), (b) data, (c) generic t + t MC events 
and (d) MC events on the vs cos 

^miss— thrust plane. Selected regions are indicated by 

the curves with arrows. 



to avoid fake 7/ candidates from 7r° — > 77. We remove events if a 7r° with invariant mass in the 
range 0.10 GeV/c 2 < M 77 < 0.16 GeV/c 2 is reconstructed by a photon from the 7/ candidate 
and another photon with E 1 > 0.05 GeV. Figure E] shows the p — > 7r + 7r~ and 77' — > p7 mass 
distributions. Dominant backgrounds for this mode come from r™ — > h~ p°v T (+7r°) with a 
photon from 7r° decay, beam background and ISR. As shown in Fig. we see no 7/ peak 
either in data or in MC since decay modes with an 7/ are very rare and are not included in 
the generic r decay model [l8j |. 

To suppress fake candidates from beam background and ISR, we require the photon energy 
to be greater than 0.25 GeV for the barrel and 0.40 GeV for the forward region. Similar to 
the 77 mode, we require the following: 5.29 GeV < <11.0 GeV, m tag < 1.777 GeV/c 2 , 
0.50 < cos Of^, (see Fig. El (a), (b) and (c)). The requirement on the number of the photon 
candidates for the signal is n SIG < 1, (see Fig. El (d)). We apply the following requirements 
Pmiss > — 5771^ — 0.2 and p miss > 2m 2 niss — 0.3 for hadronic tags, and p m i SS > —Sm^^ — 0.2 
and jJmiss > 1.2m^ m — 0.3 for leptonic tags, respectively. Following all the selection criteria, 
the signal detection efficiencies for the electron and muon modes are 5.28% and 6.00%, 
respectively. 
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10 

(GeV/c 2 ) 2 

FIG. 2: Scatter-plots of missing momentum, p m i ss , vs missing mass squared, Weiss' f° r ( a ) leptonic 
and (b) hadronic tags. Selected regions are indicated by lines. The data are the dots. The filled 
boxes show the signal MC (r~ — > fJ>~rj) distribution with arbitrary normalization. 



Event selection for the rj' — ► 7/(— > 77)7r + 7r mode 

For 77 meson reconstruction, we apply the same selection criteria as in the r — > £r](— > 77) 
analysis. The 77 — > 77 mass window is chosen to be 0.515 GeV/c 2 < m 77 < 0.570 GeV/c 2 and 
we reject photons from tt° decay. Next, we reconstruct 7/ candidates using an 77 candidate 
and two oppositely charged tracks consistent with being pions. We require P(e) < 0.1 
for both tracks in the 77' candidate. The 77' mass window is chosen to be 0.920 GeV/c 2 
< m m +^- < 0.980 GeV/c 2 , which corresponds to ±3. Oct. 

We apply the same cuts on cos0y^ ust _ miss , E™, cos6™ the invariant mass on the tag 
side, the number of photons on the signal side and m 2 ^ vs. p m i SS cut as on the t]' — > p7 
analysis. We also impose the following requirements: p m i SS > — 4m 2 niss — 0.8 and p m i SS > 
2.5m^ iss — 0.2 for hadronic tags, and p miss > — 3m^ iss and p miss > 1.5m^ ias — 0.5 for leptonic 
tags, respectively, where p tag is in GeV/c. Following all the selection criteria, the signal 
detection efficiencies for the electron and muon modes are 4.75% and 5.47%, respectively. 



Event selection for the 7r°(— > 77) modes 

The 7T° candidate is required to satisfy the condition 0.115 GeV/c 2 < M 77 < 0.152 
GeV/c 2 with p n o > 0.1 GeV/c on the signal side. We require the same cuts as for the 
v{~ > 77) mode since the final state is the same. We change the photon energy thresholds 
for the 7T° candidate > 0.90 GeV, E^ > 0.20 GeV and pi > 1.5 GeV/c compared to 
t — > £f](— > 77) (shown in Fig. 0(a), (b) and (c)). Furthermore, the opening angle between 
£ and n° in the CM system (cos 9f J, 1 ) should be in the range 0.5 < cosOfJ, 1 < 0.8 (shown 
in Fig. 0(d)). Following all the selection criteria, the signal detection efficiencies for the 
electron and muon modes are 4.35% and 5.03%, respectively. 



8 



I |t — >37t) 




0.5 0.55 0.6 

M 3 „ (GeV/c 2 ) 



FIG. 3: Invariant mass distribution of r/ — > 7r + 7r~7r candidates. While the signal MC 
(r~ — > fi~rj(-^>- 7T + 7r~7r )) distribution is normalized arbitrarily, the data and background MC 
are normalized to the same luminosity. The selected region is indicated by arrows between the 
vertical lines 
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FIG. 4: Kinematic distributions used in the event selection after the initial r] — > 7r + 7r~7r selection, 
muon identification and missing particles requirements: (a) The total visible energy in the CM 
frame; (b) the opening angle between the muon and rj candidate in the CM frame; while the signal 
MC (t~ — > fi~r](— > 7r + 7r~7r )) distribution is normalized arbitrarily, the data and background MC 
are normalized to the same luminosity. Selected regions are indicated by the arrows between the 
vertical lines. 



SIGNAL REGION AND BACKGROUND ESTIMATION 

Signal candidates are examined in the two-dimensional plots of the £~M° invariant mass, 
Mj nv , and the difference of their energy from the beam energy in the CM system, AE. A 
signal event should have M inv close to the r-lepton mass and AE close to zero. For all modes, 
the Mi nv and AE resolutions are parameterized from the MC distributions with asymmetric 
Gaussian shapes to account for initial state radiation and ECL energy leakage for photons. 
The resolutions in M inv and AE are given in Table 
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M„„ (GeV/c 2 ) M„ (GeV/c) 

FIG. 5: The p — ► 7r + 7r _ (left) and 77' — ► ^7 (right) mass distributions. While the signal MC 
(t~ — > fJ>~r)') distribution is normalized arbitrarily, the data and background MC are normalized 
to the same luminosity. Selected regions are indicated by the arrows. 



Mode 


(MeV/c 2 ) 


</L (MeV/c 2 ) 


off (MeV) 


(MeV) 


77) 


14.7 


19.4 


30.3 


61.4 


fJLrj(— ► 7r + 7r~7r°) 


7.2 


8.5 


18.5 


36.4 


«jH 77) 


14.0 


19.8 


37.3 


62.4 


er](— » 7T + 7r~7r°) 


7.6 


9.3 


19.4 


41.8 




7.8 


9.0 


16.8 


34.1 


fJ.rf(— ► 777r + 7r~) 


11.2 


19.1 


27.1 


53.5 


er/(-» /J7) 


9.2 


10.4 


19.6 


40.0 


er/(— > i]Tr + ir~) 


10.3 


21.9 


26.1 


59.4 




14.9 


19.1 


33.8 


63.0 


evr 


12.7 


23.1 


35.6 


64.6 



TABLE I: Summary of M inv (MeV/c 2 ) and AE resolutions (MeV) 



To evaluate the branching fraction, we use elliptical regions, which contain 90% of the MC 
signal that satisfies all cuts. We find that an elliptical signal region gives better sensitivity 
than a rectangular one. The signal regions are shown in Fig. |H1 the corresponding signal 
efficiencies are given in Table HT1 

We blind the signal region so as not to bias our choice of selection criteria. Figures El 
and EH show scatter-plots for data and signal MC samples distributed over ±10<r in the 
Mi nv — AE plane. As there are few remaining MC background events in the signal ellipse, 
we estimate the background contribution using the M- mv sideband regions. Extrapolation to 
the signal region assumes that the background distribution is flat along the M inv axis. We 
then estimate the expected number of the background events in the signal region for each 
mode using the number of data events observed in the sideband region inside the horizontal 
lines but excluding the signal region as shown in Fig. |S] and El The numbers of background 
events in the 90% elliptical signal region are also shown in Table ITT1 
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FIG. 6: Kinematic distributions used in the event selection after rj' — > ^7 selection, muon identifi- 
cation and missing particles requirements: (a) the total energy in the CM frame; (a) the invariant 
mass on the tag side; (c) the opening angle between the missing particle and tag-side track in the 
CM frame; (d) the number of photons on the signal side. While the signal MC (r~ — > [j,~r]'(—> pj)) 
distribution is normalized arbitrarily, the data and background MC are normalized to the same 
luminosity. Selected regions are indicated by arrows from the marked cut boundaries. 



Systematic uncertainties for M° reconstruction are 3.0%, 4.0%, 4.0%, 5.0% and 3.0% for 
rj — > 77, 7] — > 7T + 7r~7r°, r]' — > p'-f, 7]' — > r]ir + TT~ and it — > 77, respectively. Furthermore, the 
uncertainties due to the branching ratios of the M° meson are 0.7%, 1.8%, 3.4% and 3.5% for 
rj — > 77, r] — > 7T + 7r~7r°, r]' — > /ry and r/ — > ?77r + 7r~, respectively 0|. For the n° veto we take 
a 5.5% uncertainty for rj — > 77 while a 2.8% uncertainty is assigned to the 77' — > p7 mode. 
The uncertainties in the trigger (0.5—1.0%), tracking (2.0%), lepton identification (2.0%), 
MC statistics (1.0—1.5%), luminosity (1.4%) are also considered. All these uncertainties are 
added in quadrature, and the total systematic uncertainties are shown in Table ITT1 

While the angular distribution of signal r decays is initially assumed to be uniform in 
this analysis, it is sensitive to the lepton flavor violating interaction structure The spin 
correlation between the r lepton on the signal and that on the tag side must be considered. 
A possible nonuniformity was taken into account by comparing the uniform case with MC's 
assuming V — A and V + A interactions, which result in the maximum possible variations. 
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FIG. 7: Kinematic distributions used in the event selection: (a) higher energy and (b) lower energy 
of a photon from the ir° candidate (Eyi and E^)', (c) momentum of a muon (p^); (d) the cosine 



of the opening angle between the muon and tt in the CM frame (cos 6^ ). While the signal MC 
(r~ — > /U - 7T°) distribution is normalized arbitrarily, the data and background MC are normalized 
to the same luminosity. Selected regions are indicated by arrows from the marked cut boundaries. 



No statistically significant difference in the M inv - AE distribution or the efficiencies is found 
compared to the case of the uniform distribution. Therefore, systematic uncertainties due 
to these effects are neglected in the upper limit evaluation. 

We open the blind and find no data events in the blinded region after event selection. 
Only in the case of r — ■> /i7r°(^ 77), one event is found in the elliptical region. Since no 
statistically significant excess of data over the expected background in the signal region 
is observed, we set upper limits for branching fractions. The upper limit on the number 
of signal events at the 90% confidence level (C.L.) sq n in cluding systematic uncertainty is 
obtained with the use of the Feldman- Cousins method |20J] calculated by the POLE program 
without conditioning j23|. The upper limit on the branching fraction {B) is then calculated 
as 
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Mode 


R „ 


e 


»o 


s 


Total Sys. 


S90 




T 






fi 49% 


n 4n-i-n 9Q 


n 

u 


7.1% 


9 1 


T 




/• + — 0\ 

/W7(— >• tt^tt 7T J 




fi 84% 


D 94-1-0 94 


n 


5.6 


% 


9 9 




T 


-> er?(^ 77) 




4 57% 


n 9^-l-n 9^ 


n 

u 


7.1% 


9 9 
_ . _ 


T 




e?y( — ^ tt^tt tt^) 


0.226 


4.72% 


53±0 53 





5.6 % 


2.0 




T 


—> \XT] (— > P7J 


VJ.ZicfO ALU 


^ 4D% 




n 

u 


6.8 


% 


9 9 


r 






0.443x0.3943 


4.92% 


— U.U 





8.9 


% 


2.5 




T 


-» er/(-> ^7) 


0.295x1.0 


4.76% 


o+ - 33 
u - u -o.o 





6.8 


% 


2.5 


T 




erf{— > r]ir + 7r~) 


0.443x0.3943 


4.27% 


o +0 - 24 





8.9 


% 


2.5 




T 


-> A»t (-» 77) 


0.98798 


4.53% 


0.58±0.34 


1 


4.5 


% 


3.8 




T 


— > e7r°(^ 77) 


0.98798 


3.93% 


0.20±0.20 





4.5 


% 


2.2 



TABLE II: Results of the final event selection for the individual modes: B^fo is the branching 
fraction for the M° decay; bo and s are the number of expected background and observed events 
in the signal region, respectively; "Total sys." means the total systematic uncertainty; Sgo% is the 
upper limit on the number of signal events including systematic uncertainties. 



Mode 


M° subdecay mode 


Upper limit of at 90%C.L. 




77 -> 77 

7/ — » 7T + 7r~7r° 

Combined 


1.2xl0~ 7 
2.0xl0~ 7 
6.5xl0~ 8 


t~ — ► e~rj 


7/ -> 77 

7/ — » 7T + 7r~7r° 

Combined 


1.7xl0" 7 
2.6xl0~ 7 
9.2xl0" 8 




7/ -> p7 

7/ — > 7/7r + 7r~ 

Combined 


1.9xl0~ 7 
4.1xl0~ 7 
1.3xl0" 7 


t~ — ► e~rj 


77 -»-p7 

7/ — » 1]TT + TT~ 

Combined 


2.5xl0~ 7 
4.7xl0~ 7 
1.6xl0~ 7 


T~ — > /J~7T° 


7T° —5- 77 


1.2xl0~ 7 


r~ — > e~7r° 


7T° — ► 77 


8.0xl0~ 8 



TABLE III: Summary of upper limits on B at 90%C.L. 



where Bm° is taken from PDG 18|| and N TT = 357.7 x 10 6 is the number of r— pairs produced 
in 401 fb -1 of data. We obtain N TT using o TT = 0.892 ± 0.002 nb, the e + e" — > t + t~ cross 
section at the T(4S") resonance calculated by KKMC [l^|. The upper limits for the branching 
fractions are in the range B{t~ — > £~M°) < (6.5 — 16) x 10~ 8 at the 90% confidence level, 
respectively. A summary of the upper limits is given in Table IIHI These results improve the 
previously published limits 0] by factors of 2.3—6.4. 
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FIG. 8: Scatter-plots of data in the M- mv - AE plane: 
±10cr area for the r~ — ► /i~r}(—> 77), r~ — > e~rj{-^ 77), r" 



1.75 



^(Gev/cf 



a), (b), (c) and (d) correspond to the 



7T 7T 7T 



and r" 



7T 7T 7T 



0> 



modes, respectively. The filled boxes show the MC signal distribution with arbitrary 
normalization. The elliptical signal region shown by the solid curve is used for evaluating the signal 
yield. The region between the horizontal lines excluding the signal region is used to estimate the 
expected background in the elliptical region. 



DISCUSSION 

The branching ratio for the r~ — > p,~7] mode is enhanced by Higgs-mediated LFV if large 
mixing between a left-hand scalar muon and scalar tau in the corresponding SUSY model 
occurs (H and can be written as 

7 /tan/3\ 6 ZlOOGeV/cV 
B(T -^^) = 8.4X10-^] [—^-) , CO 

where is the pseudoscalar Higgs mass and tan f3 is the ratio of the vacuum expectation 
values. From our upper limit on the branching fraction for the t~ — > pTr\ decay, some region 
of mA and tan j3 parameters can be excluded. Figure [TO] shows the excluded region in the 
rriA vs tan/3 plane. Figure ITU1 also shows the constraints at 95% C.L. from the CDF 22l |. 
D0 @] and LEP2 experiments Q. The excluded regions from the CDF, D0 and LEP2 
experiments are shown with fi > and m™ ax . Our result has a sensitivity competitive with 
that of the CDF and D0 experiments, which searched for pp — > (pb(b) and t + t~ events, 
where is a neutral Higgs boson in the minimal supersymmetric standard model (0 = h,H 
and A). 
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FIG. 9: Scatter-plots of data in the M\ m - AE plane: (e), (f), (g), (h), (i) and (j) correspond to the 
±10<7 area for the r~ — > /i - r/(— > pj), r~ — > e~r/(— > pj), t~ — > n~r]'(—> 7]ir + ir~), r~ — > e _ r/(— > 
t]7t + tt~), r~ — > /i~7r°(^ 77) and r~ — > e~7r°(— > 77) modes, respectively. The data are indicated 
by the solid circles. The filled boxes show the MC signal distribution with arbitrary normalization. 
The elliptical signal region shown by the solid curve is used for evaluating the signal yield. The 
region between the horizontal lines excluding the signal region is used to estimate the expected 
background in the elliptical region. 



SUMMARY 

We have searched for lepton flavor violating r decays with a pseudoscalar meson (rj, rj' 
and 7T°) using 401 fb _1 of data. No signal is found and we set the following upper limits 
of the branching fractions: B(t~ — > e~rj) < 9.2 x 10~ 8 , B(t~ — > ^ff) < 6.5 x 10~ 8 , 
B(t~ -> e~rf) < 1.6 x 1(T 7 , B(t~ -> n'lf) < 1.3 x 10~ 7 B{r~ -> e-vr ) < 8.0 x 10~ 8 
and B(t~ — > /U _ 7r°) < 1.2 x 10~ 7 at the 90% confidence level, respectively. These results 
improve the previously published limits by factors from 2.3 to 6.4 and help to constrain 
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FIG. 10: The excluded region in the rriA vs tan/3 plane from our results at 90% C.L. and other 
experiments at 95% C.L. from CDF Q, D0 Q, LEP 0. The excluded regions from the CDF, 
D0 and LEP2 experiments are shown with [i > and m™ ax . 



physics beyond the Standard Model. 
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